Nematic, vector-multipole, and plateau-liquid states in the classical 0(3) pyrochlore 
antiferromagnet with biquadratic interactions in applied magnetic field 
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The classical bilinear-biquadratic nearest-neighbor Heisenberg antiferromagnet on the pyrochlore lattice 
does not exhibit conventional Neel-type magnetic order at any temperature or magnetic field. Instead spin 
correlations decay algebraically over length scales r < ~ VT, behavior characteristic of a Coulomb phase 
arising from a strong local constraint. Despite this, its thermodynamic properties remain largely unchanged 
if Neel order is restored by the addition of a degeneracy-lifting perturbation, e.g., further neighbor interac- 
tions. Here we show how these apparent contradictions can be resolved by a proper understanding of way in 
which long-range Neel order emerges out of well-formed local correlations, and identify nematic and vector- 
multipole orders hidden in the different Coulomb phases of the model. So far as experiment is concerned, our 
results suggest that where long range interactions are unimportant, the magnetic properties of Cr spinels which 
exhibit half-magnetization plateaux may be largely independent of the type of magnetic order present. 
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PACS numbers: 75.10.-b, 75.10.Hk 75.80.+q 



I. INTRODUCTION 

Frustrated magnets have long been studied as a paradigm 
for complex behavior in condensed matter and statistical 
physics^. The most widely studied systems are frustrated an- 
tiferromagnets (AF), where competing interactions suppress 
classical Neel order In some highly frustrated magnets, spins 
do not order at any temperature, and the ground state retains 
only very short-ranged spin-spin correlations. The resulting 
state is generally termed a "spin liquid". 

However it is also possible to invert this paradigm and think 
of highly frustrated magnets as systems where local "order" is 
robust enough to survive, even where long range order has 
been obliterated by fluctuations. Conventional Neel order can 
then easily be restored — albeit with a relatively low critical 
temperature — by any perturbation which forces long-range 
coherence on this preformed local order Moreover, where 
quantities such as heat capacity and magnetic susceptibility 
are controlled by local fluctuations, the thermodynamic prop- 
erties of the globally-disordered spin-liquid phase may be 
practically indistinguishable from those of the magnetically 
ordered phase. 

In this paper we explore this "bottom-up" formulation of 
frustration, showing how the different multipolar and spin liq- 
uid states of a simple classical frustrated antiferromagnet in 
applied magnetic field already contain the seeds of long range 
Neel order. The model which we consider is the antiferro- 
magnetic nearest-neighbor Heisenberg model with additional 
biquadratic interactions b 

n = [s,.Sj-6(s,.Sj)2] -h-^s,, (1) 

where the sum {ij) i runs over the nearest neighbor bonds of a 
pyrochlore lattice (Fig.[T]). All the energy scales including h = 
|h| and temperature T are measured in units of Ji hereafter 
This model was introduced in Ref 101 to explain the dra- 
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FIG. 1: (Color online) 16-site cubic unit cell of the pyrochlore lattice 
— a network of corner-sharing tetrahedra. Exchange interactions Ji 
are associated with the first neighbor bonds {ij}i, and J3 with the 
two (inequivalent) types of third-neighbor bond 



matic half-magnetization plateau observed in Cr spinels^—. 
In this case the biquadratic interaction b originates in a strong 
coupling to the lattice. However such terms can also be of 
electronic origin, and quite generally they can be taken to 
characterize the effects of quantum and/or thermal fluctua- 
tions in a frustrated magnet^Ji. Thus we anticipate many of 
our results will also be relevant for the quantum model. Re- 
cent results for the 5 = 3/2 XX Z pyrochlore AF in applied 
magnetic field suggest that this is indeed the cas o'^''^ . 

It is well known that the classical Heisenberg model with 
the nearest-neighbor bilinear couplings only does not exhibit 
Neel-type magnetic order on the pyrochlore lattice at any 
temperature'-'- . As we shall see, these arguments are es- 
sentially unchanged by the introduction of magnetic field, or 
by nearest-neighbor biquadratic interactions b. The system 
can however be brought to order by introducing an interaction 
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FIG. 2: (Color online) (a) Magnetic phase diagram of the classical 
pyrochlore antiferromagnet with biquadratic interactions b — 0.1, 
and additional third-neighbor interactions J3 — —0.05, as deter- 
mined by classical Monte Carlo simulatior*— The form of four- 
sublattice Neel order is illustrated, together with the irreducible rep- 
resentation (irrep) of the tetrahedral symmetry group Td to which it 
belongs. For these parameters, the model provides a good description 
of the half-magnetization plateau seen in CdCr20^^-. (b) Equivalent 
magnetic phase diagram in the absence of any longer range interac- 
tions. Canted Neel states are replaced by phases with multipolar or- 
der, while the collinear half-magnetization plateau state gives way to 
a collinear spin liquid. Crosses denote the crossover at temperature 
T* from paramagnet to the plateau-liquid state, as determined by a 
peak in the heat capacity. In both (a) and (b) circles with solid (red) 
lines denote first-order phase boundaries, while those with dashed 
(blue) lines denote second-order ones. Both T and h are measured 
in units of Ji . 



which links spins in different tetrahedra, for example, 

T^LRO ^ J3 ^ s, • S, , (2) 

where {ij)3 runs over the two (inequivalent) sets of third 
neighbor bond shown in Fig.[T] 

For ferromagnetic (FM) J3 < 0, this specific form of Ti}-^'^ 



leads to the four-sublattice long-range order (LRO) described 
in Ref. 10], and to the finite temperature transitions shown in 
Fig. [2ta) '^''^ . Four sublattice order can also be stabilized by 
AF second neighbor interaction J2 . More generally, however, 
the type of order which results depends on the details of the in- 
teraction 7^LRqi8ji9 -j^g system can therefore be tuned at will 
between different types of ordered state, simply by changing 
n^^'^. From this we conclude that, as a function of magnetic 
field h, for 'H^^'^ = 0, there must be a line of second-order 
multicritical — or first-order multifurcative points — separat- 
ing a huge set of different ordered phases. 

The main purpose of this paper is to explore the symmetry 
breaking which persists in the limit of Ti}-^'-' — > for finite 
biquadratic interaction h and finite temperature T. In order to 
make the problem accessible to large scale Monte Carlo (MC) 
simulation, we consider the classical 5 = |S| ^ 00 limit of 
Eq. ([T]i, reseating variables such that S =1. 

Using a mixture of classical MC simulation, analytic low- 
T expansion, and simple field theoretical arguments, we find 
a set of phases in the h-T plane which exhibit power-law de- 
cay of spin correlation functions. Two of these phases possess 
long-range nematic or vector-multipole order and, most inter- 
estingly, the magnetization plateau persists in the absence of 
conventional magnetic order. We show how all of these results 
can be understood — and even anticipated — from a proper 
understanding of the geometry of the pyrochlore lattice, and 
the way in which a single tetrahedron behaves in magnetic 
field. Our findings are summarized by the h-T phase diagram 
shown Fig. |2tb). 

So far as experiment is concerned, our main conclusion will 
be that the thermodynamic properties of the pyrochlore anti- 
ferromagnet in applied magnetic field are mostly determined 
by symmetry breaking at the level of single tetrahedron. Local 
order is well-formed for 'H}-^'~' = 0, and many properties of 
the system are therefore insensitive to the details of the LRO 
order present. Thus the very simple phase diagram derived in 
Ref. J^l and its finite temperature generalization in Refs. ifTill 
and 1I15I1 [reproduced in Fig.|2a)], are applicable for a wide 
variety of different 'H}-^'~' . 

The paper is structured as follows: In Sec. |ll] we briefly 
review the basic physics of the Heisenberg model on the py- 
rochlore lattice. Definitions are given of order parameters for 
conventional Neel (dipolar) order, and of rank-two tensor or- 
der parameters which can be used to signal multipolar order 

Then, in Sec. |lll] we use these tools to construct the 
h-T phase diagram of the pyrochlore AF with additional bi- 
quadratic interactions [Fig. |2jb)]. Thermal fluctuations pre- 
serve the extensive degeneracies present in the ground state, 
and fail to select any conventional long-range dipolar or- 
der Despite this, the thermodynamic properties of the sys- 
tem and the topology of the phase diagram are essentially un- 
changed — the magnetization plateau survives and nematic 
and vector-multipole phases corresponding to the two differ- 
ent canted states are shown to exist for fields below and above 
the magnetization plateau [Fig.|2a)]. 

In Sec. |IV] we explore the way in which long-range Neel 
order is recovered as a FM third-neighbor interaction J3 is 
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FIG. 3: (Color online) (Left panel) Gound-state phase diagram of a single classical tetrahedron as a function of magnetic field h and di- 
mensionless coupling constant b, taken from Ref. 0]. Solid (red) lines denote first and dashed (blue) lines second order transitions. Spin 
configurations and relevant irreducible representations (irreps) are shown in each case. (Right panel) Symmetries of the Ai, E, and T2 irreps 
of the tetrahedral group Td used to classify different states. Solid (red) lines have negative weight; hollow (blue) lines have positive weight. 
Thin (black) lines have zero weight. See also Eq. Q. 



"turned on", focusing on the half-magnetization plateau for 
/i w 4. For small | J3I, the system now exhibits two charac- 
teristic temperature scales — an upper temperature T* k, b 
at which the gap protecting the magnetization plateau opens, 
and a lower temperature Tn « (!)(| J3I) at which the system 
exhibits long-range magnetic order. This is contrasted with 
the situation for h ^ Q, where the system also exhibits two 
characteristic temperature scales, but these correspond to suc- 
cessive phase transitions : a nematic transition at Tq ^ b and 
a Neel ordering at Tn ^ 0(| J3I). We discuss the nature of 
these transitions for J3 — > 0, identifying a line of first-order 
multifurcative points at J3 = 0. And, for J3 = 0, we identify 
an unusual continuous transition from the coulombic plateau 
liquid to the vector-multipole phase. At low temperatures this 
transition appears to have mean field character. 

Finally, in Sec. |V] we conclude with a discussion of the 
broader implications of these results. 



II. DEGENERACIES IN FINITE MAGNETIC FIELD 

A. Geometrical arguments 

The pyrochlore lattice (Fig.[T]i is the simplest example of a 
three-dimensional (3D) network of corner sharing complete 
graphs. Its elementary building block is the tetrahedron, in 
which every site is connected to every other site, i.e. the tetra- 
hedron is a complete graph of order four. Tetrahedra in the py- 
rochlore lattice can be divided into A and B sublattices, with 
each lattice site shared between an A- and a B-sublattice tetra- 
hedron. The centres of the two types of tetrahedra together 
form a (bipartite) diamond lattice The overall symmetry of 
the lattice is cubic. 



As such, the pyrochlore lattice is a natural 3D analogue of 
the 2D kagome lattice, a corner sharing network of triangles 
(complete graphs of order three). In fact the [111] planes 
of the pyrochlore lattice are alternate kagome and triangu- 
lar lattices, composed of the triangular "bases" of tetrahedra 
and their "points", respectively. Much of the unusual physics 
of the kagome lattice also extends to its higher dimensional 
cousin. 

Lattices composed of complete graphs have the special 
property that bilinear quantities on nearest neighbor bonds can 
be recast as a sum of squares. Thus for 6 = the Hamiltonian 
Eq. ([T]i can be written 

H = 4^(m-^) -^+const., (3) 

tetra ^ ' 

where the sum runs over tetrahedra, and 

M= i(Si+S2 + S3 + S4) (4) 

is the magnetization (per site) of a given tetrahedron. For 
/i = 0, a simple classical counting argument shows that two of 
the eight angles needed to determine the orientation of the four 
spins in any given tetrahedron remain undetermined. Nearest 
neighbor interactions do not select one unique ground state on 
the pyrochlore lattice but rather the entire manifold of states 
for which |M| = in each tetrahedron. Thus at T = 0, the 
system is disordered. For fields h < h^at = 8 this conclu- 
sion is unaltered by the presence of magnetic field. In this 
case the manifold of ground states is determined by the con- 
dition M = h/8 in each tetrahedron, and the magnetization 
is linear in h up to the saturation field /igat = 8. (We recall 
that magnetic field is measured in units of Ji, so that in fact 

hsat — 8J1.) 
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In order to understand how nearest-neighbor biquadratic 
interactions b select among this manifold of states, it is suf- 
ficient to solve the problem of a single tetrahedron embed- 
ded in the 3D lattice. This problem was considered in 
Ref For 6 > 0, biquadratic interactions select coplanar 
(and coUinear) configurations from the larger ground-state 
manifold of Eq. There are four dominant phases, illus- 
trated in Fig.|3]: 

(i) a 2:2 coplanar canted state for low field 

(ii) a 3:1 collinear (uuud) half-magnetization plateau state 
for intermediate field 

(iii) a 3:1 coplanar canted state for fields approaching satu- 
ration 

(iv) a saturated (uuuu) state for large magnetic field 

h > /isat 

An exhaustive enumeration of possible states is given in 
Ref 

Up to this point, we have not been specific about how the 
tetrahedron was embedded in the lattice. It could, trivially, 
form part of a state with Neel order, e.g., the simple four- 
sublattice order favored by FM J3. However, there are in- 
finitely many other ways of joining 2:2 or 3:1 tetrahedra to- 
gether at the corners, and not all of them correspond to Neel 
ordered states. In fact the ground state manifold retains an ex- 
tensive Ising-like degeneracy for all h < hgat, and as a result 



the system remains "disordered". The nature of this degener- 
acy, and its consequences, are explored in some detail below. 



B. Bond order parameters 

Where Neel order is present, it can be detected in the re- 
duced spin-spin correlation function 

Dir,,) = (S, • Sj) - . (5) 

Here is the expectation value of the squared magnetization 
per spin, 

i 

which vanishes in the absence of magnetic field. N is the total 
number of spins. The simplest form of order supported by the 
pyrochlore lattice is the four-sublattice Neel order favored by 
FM J3, as illustrated in Fig.|2a). 

Written in terms of the minimal four-site unit cell of the py- 
rochlore lattice, four-sublattice order has momentum q = 0, 
and different states can easily be classified using the Ai, E, 
and T2 irreducible representations (irreps) of the symmetry 
group Td for a single tetrahedron : 
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(7) 



where the spins belong to a single tetrahedro n"!^"'^' — cf. 
right panel of Fig. |3] We can use these irreps to define bond 
order parameters 

tetra 

and associated generalized susceptibilities 

Xf^' = f [((Ar^?)-(Af^=')^], (9) 

where the sum X^tetra ^'^^ '^^^^ ^ 1^ independent A- 
sublattice tetrahedra, and A^, is the vector associated with 
the V = {E, T2} irreps of the tetrahedral symmetry group Td, 
namely, Ae = (AEa,AE,2) and Aj^ = (Aj^.i, At2,2, Aj^.s)- 
These order parameters allow us to distinguish the sharp 
first-order transition between orders in the T2 and E irreps. 



and the transition from the Ai to T2 states at high field — cf. 
Fig. 12 a) — but not the more subtle second-order transition 
between the T2 symmetry uuud plateau state and the T2 sym- 
metry 3:1 canted state. These are none the less distinct phases 
— the collinear and canted T2 states are connected by a zone- 
center (i.e., q = 0) excitation which is gapped in the collinear 
uuud state, and becomes soft at the critical field marking the 
onset of the 3:1 canted stata^S. The condensation of this soft 
spin mode corresponds to the emergence of order in the trans- 
verse spin components (5^), {S"^) — i.e., canting of spins 
away from the z axis (the direction of applied magnetic field). 

The bond order parameters defined in Eq. (|7]i couple di- 
rectly to the lattice, through changes in bond length^. They are 
therefore particularly well suited to describing simple Neel or- 
dered states, where the magnetic ordering is driven by the lat- 
tice effects. However the irreps on which they are based also 
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provide a useful measure of the correlation which survives in 
the absence of long range order, a central question for this pa- 
per To this end, we introduce a measure of local correlation 



^local 



tetra 



and its associated generalized susceptibility 



xl°"' = ^[((AL°"'n-(AL°=^'> 



l\21 



(10) 



(11) 



For a single tetrahedron, X^°"^ and A^'"'^^' are identical. On 
a lattice, A^'^^' lacks crossterms between different tetrahedra 
present in Xf°'°^\ and is therefore a measure of correlation in 
the absence of long range order. We return to these points 
below. 



C. Rank-two tensor order parameters 

Not all of the phases supported by the Hamiltonian Eq. ([T]i 
can be described using the bond order parameters Eq. (|7]i. In 
Appendix lAl we formally classify the different types of sym- 
metry breaking which can arise in this model at the level of a 
single site. Here we restrict ourselves to the simplest possible 
generalization from Neel to multipolar order; both the T2 and 
E symmetry canted states possess order of transverse (i.e., x 
and y) spin components which vanishes in the collinear uuud 
state, and which can survive even in the absence of conven- 
tional (canted) Neel order. 

To describe this, it is convenient to introduce the rank-two 
tensor order parameters 
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where the local quadrupole moments 
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are summed over all lattice sites i. 

Where spin rotational symmetry is not akeady broken by 
magnetic field, i.e., for /i = 0, spins may select a common 
axis without selecting a direction on it. This is conventional 
nematic order, of the type exhibited by uniaxial molecules, 
and can be detected using the order parameter 

+ (g"^)' + (Q"^)' + (0'^")' (18) 

which is invariant under 0(3) rotations. This order param- 
eter takes on its maximal value (Q^) 4/3 in a perfectly 
collinear state, such as the 2:2 state for T 0. 



TABLE I: Classification of tensor order operators according to ro- 
tational symmetry about a z axis defined by magnetic field: Each 
forms an irrep transforming like e*"'*, where n is an integer and (j) is 
the polar angle in the xy plane. Also indicated are the finite values 
of the order parameters in the 2:2 and 3: 1 canted states. 



order par. 



tensor operators 



2:2 



1:3 



finite 



finite 



finite 




finite 
finite 



finite 
finite 



In what follows we will also make use of the correlation 
function measuring collinearity 



(19) 



considered in Ref. IH. As defined, -1/2 < (P(ry )) < 1, 
taking on the value {Pijij)) = for uncorrected spins. In 
fact P{vij) can also be expressed in terms of quadrupolar op- 
erators as 



and it follows that 



iV2 



(20) 



(21) 



At finite h, the 0(3) invariant correlation function Eq. (fT9b 
still provides a useful measure of collinearity, but does not by 
itself signal a broken symmetry. In this case it is convenient 
to group quadrupoles according to way in which they trans- 
form under the remaining 0(2) rotations about the direction 
of magnetic field — conventionally the z axis. We therefore 
consider 



(22) 

(23) 
(24) 



where the magnetic field is assumed to be parallel to the 
z axis. Each of the separate irreps Q-'-'" transforms like 
{cos 770, sin 770} — or equivalently, e™"^ — where (j) is the 
polar angle in the plane perpendicular to the magnetic field. 
They can therefore be used as order parameters to detect the 
77-fold breaking of rotational symmetry in the xy plane. The 
conventional nematic order parameter with full 0(3) symme- 
try, Eq. (fTSl l, is given by the sum of squares 



IQ' 



IQ' 



(25) 



In finite magnetic field, the one-dimensional irrep Q^'^ 
does not contain any information about broken symme- 
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tries and can generally be discarded. However the two- 
dimensional irreps Q-'-^^ and Q-'-^^ distinguish different or- 
dered phases. In the 2:2 canted phase the mean square value 
of Q-'--^ takes on a finite value 

(|Q^.2|2^ ^ ^(g.^-,^)2^(QX.,)2^ ^26) 

This is another form of nematic order of the transverse spin 
moments — one transforming like e*^*^ — and reflects the fact 
that spins select a common plane in which to cant. At the same 
time mean square value of Q-'-'^ — which transforms as e*"^, 
i.e., a vector in the xy plane — vanishes. 

Similarly, (IQ^'^P) takes on a finite value in the 3:1 canted 
phase. However in this case the 3:1 asymmetry of the canted 
spin configuration defines a direction in the xy plane, and 

(IQ'-'T) = ((g"")' + (Q^^)'> > (27) 

is also finite. The 3: 1 canted phase therefore possess a form of 
vector-multipole order. These facts are summarized in Table|I] 
In what follows we concentrate almost exclusively on 
phases which do not exhibit conventional magnetic order, as 
defined by D{r) in Eq. (|5]l, and characterize these states using 
the rank-two tensor order parameters listed in Table J] For 
further details of conventional Neel phases, and comparison 
with experiment, we refer the interested reader to Ref. ifTsll . 
Rank-three tensors which also occur as order parameters in 
the present model are discussed in Appendix iB] 



D. General considerations 

Many frustrated systems with disordered ground states 
manage none the less to order at finite temperature. This effect 
is known as "order from disorder" and occurs where there is a 
net entropy gain in selecting one particular state out of the dis- 
ordered manifold. Entropy is gained where a given spin con- 
figuration (typically, collinear or coplanar) has a higher den- 
sity of low-energy excitations than its peers. However this 
entropy gain must be sufficient to offset the entropy lost by 
choosing one state out of the manifold. Where the ground 
state manifold has an extensive degeneracy, this is a very 
strong constraint. Order-from-disorder effects are known to 
select one particular Neel ordered ground state in e.g., the 
frustrated square lattice^, but fail to do so in the case of the 
more frustrated kagome lattice^^''^'*. 

Even where fluctuations fail to stabilize one particular Neel 
ground state, they can still select a subset of states from the 
ground state manifold with a smaller — but none the less ex- 
tensive — degeneracy. This subset (submanifold) of states 
will not exhibit the long range spin-spin correlations which 
are the hallmark of conventional Neel-type magnetic order 
However this does not necessarily mean that the system is 
truly disordered — it may well exhibit long range order of 
a more complex type. 

A good example of this second type of order-from-disorder 
effect is provided by the nearest-neighbor classical XY 
model on the pyrochlore lattice, where thermal fluctuations 



lead to nematic order with broken spin-rotational symmetry, 
but power-law decay of spin-spin correlationsii. 

In what follows we use the order parameters defined in 
Sec. Ill CI to identify phases of Eq. ([T]i which exhibit nematic 
order in the absence of Neel order We focus chiefly on dif- 
ferent forms of unconventional order found in magnetic field. 
Closely related studies in magnetic field have been made of 
the classical Heisenberg model on a kagome lattices^, and 
classical XY model on a checkerboard latticeiS.. In both these 
cases unconventional order is stabilized by thermal fluctua- 
tions. Another type of unconventional order for the pyrochlore 
lattice with FM second-neighbor interactions J2 < and 
h = was recently studied in Ref. ll26l] . In our case the main 
driving force is not fluctuations but finite biquadratic interac- 
tion b; results for order stabilized by thermal fluctuations at 
finite h and J3 but & = will be presented elsewhera^i. 

III. PARTIAL LIFTING OF DEGENERACY IN FINITE 
MAGNETIC FIELD 

A. Collinear nematic phase for h — 

In the absence of magnetic field, the ground state of Eq. ([T]) 
is determined by the conditions that (i) the total magnetization 
of each tetrahedron be zero, to minimize the antiferromagnetic 
exchange interaction Ji, and (ii) all spins be collinear, to min- 
imize the biquadratic interaction b. These conditions select an 
extensive manifold of 

flo « 1.5^/2 « 1.22^ (28) 

states with exactly two-"up" and two-"down" spins (uudd) 
in each tetrahedron. The degeneracy of this ground state man- 
ifold is of the same form as that encountered in Pauling's the- 
ory of water ice--, and we therefore refer to it as the "ice" 
manifold below. Since each spin is shared by two neighbor- 
ing tetrahedra, "up" and "down" spins form unbroken loops 
as shown in Fig.lD We return to this point below. 

The fact that the direction along which "up" and "down" 
spins point is not determined by the Hamiltonian implies that 
spin rotational symmetry must be broken spontaneously (for 
simplicity, we none-the-less to use "up" and "'down" to de- 
note the oppositely oriented spins). This can be seen in the 
spin collinearity Eq. ( fT9] l. which takes on the maximal value 
P(r) = 1 for all states in the ice manifold, implying that 
the ground state manifold has nematic (i.e., quadrupolar) or- 
der (This is explicitly confirmed by MC simulations below.) 
However, as already stated, the ground state manifold does not 
possess Neel order of any form. 

In fact it is possible to calculate the asymptotic form 
of spin-spin correlations in the ice manifold by mapping 
them onto configurations of a notional electric (or magnetic) 
field22ri^. The condition that every tetrahedron has exactly 
two-"up" and exactly two-"down" spins translates into a zero 
divergence condition for the electric (magnetic) field, and 
spin-spin correlations take on a dipolar form 
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FIG. 4: (Color online) (a) Illustration of coUinear E-symmetry ne- 
matic state at h — 0, showing loop-like coordination of parallel spins 
associated with the "ice" manifold, (b) Four-sublattice long range 
order with E symmetry induced by FM J3. (c) Canted nematic state 
with partial magnetization under applied field. 



dictated by this effective electrodynamics. This power-law 
decay of spin correlations is a signal property of the "ice" 
manifold. However it should not be taken to imply that "up" 
and "down" spins are entirely uncorrected. Within each uudd 
tetrahedron, each spin has twice as many AF aligned neigh- 
bors as FM aligned ones, and the net correlation on nearest 
neighbor bonds is 



(30) 




0.16 



FIG. 5: (Color online) Spin collinearity P [Eq. l |I9l l1 at ft = for 
b = 0.1, showing the onset of nematic order at Tq « 0.13. P is 
measured for the farthest spin pair along the (110) chains in the py- 
rochlore lattice in each system size ranging from L = 4 to L = 16, 
and averaged over the (110) chains running in different directions. 



of the tetrahedral symmetry group Td, defined in Sec. IIIBI 
and that local fluctuations of order A^^^' take on their maximal 
value A^°=^' = 16/3. 

This concludes our survey of symmetry breaking for T ~ 0, 
but it leaves open the question, what happens at finite tem- 
perature? By analogy with ordered systems where "order 
from disorder" is effective, thermal fluctuations might be ex- 
pected to select a single configuration from the "ice" mani- 
fold, and so restore Neel order To address this question, we 
have performed extensive Monte Carlo simulations using a 
local-update Metropolis algorithm to sample spin configura- 
tions. We typically perform 10^ MC samplings for measure- 
ments after 10^ steps for thermalization. We have checked the 
convergence of the results by comparing those for different 
initial spin configurations. In particular, to minimize the hys- 
teresis associated with first order transitions, we used mixed 
initial conditions in which different parts of the system are 
assigned different ordered or disordered states^. Where the 
acceptance rate in MC updates becomes extremely slow, we 
used the exchange MC method^!, to avoid local spin-freezing 
at low temperatures. Results are divided into five bins to es- 
timate statistical errors by variance of average values in the 
bins. The system sizes in the present work are up to i = 16, 
where L is the linear dimension of the system measured in 
the cubic units shown in Fig.lT] i.e., the total number of spins 
N is given by 16L^. We show the results for b = 0.1 and 
b = 0.6, case by case, both of which exhibit qualitatively the 
same behavior 

In the absence of applied magnetic field, the spin collinear- 
ity P{r) grows sharply below a transition temperature 
Tq Ri b, as illustrated in Fig. |5] As expected, for T — > 0, 
P(r) 1, implying that all spins have a single common 
axis. The nematic order parameter Q [Eq. ( fTsl lI is plotted 
in Fig. |6la), together with the heat capacity in Fig. |6tb) for a 
range of system sizes from L = 4 to L = 16. Here the heat 
capacity is calculated by the fluctuation of internal energy as 



Locally, order is well formed. More formally, we can state that 
these uudd tetrahedra belong to the two-dimensional E irrep 
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FIG. 6: (Color online) (a) Nematic order parameter Q [Eq. 1II8H. 
showing the onset of nematic order at Tq ~ 0.13; (b) Heat capacity; 
(c) Absence of long range four-sublattice order ^l'"*"^' is accompa- 
nied by well-formed local correlations X^^"'; (d) The associated lo- 
cal susceptibility shows a sharp jump at Tq, where spins in tetrahe- 
dra with preformed local order gain energy by selecting long range 
collinearity. All data are for ft = 0, 6 = 0.1, and for system size 
ranging from L = 4toL = 16. 



The sharp onset of order and jump in heat capacity imply a 
first order phase transition at Tq w 0.13. 

Treating this nematic order at the level of a Ginzburg- 
Landau theory, the free-energy terms allowed by lattice and 
spin rotational symmetries are 



J" = aQ^ + bQ^+cQ^ 



(32) 




FIG. 7: (Color online) (a) Spin correlations D [Eq. (|5}] at h — 
for b — 0.6, for r < L measured along (110) chains in the py- 
rochlore lattice and averaged over different chain directions as in 
Fig. [5] r is in units of the distance between nearest neighbor spins, 
(b) Characteristic 1 /r'' power-law decay of spin correlations in the 
low-temperature nematic phase for T < Tq ~ 0.5 associated with 
the ice manifold of 2:2 states, plotted on a log-log scale, (c) Expo- 
nential decay of spin correlations in the high-temperature paramag- 
netic phase for T > Tq ~ 0.5, plotted on a log-linear scale. In (b) 
and (c), grey lines are guides for the eye, ~ 1/r'' and ~ exp(— r), 
respectively. In all figures the data are for temperatures ranging from 
T = 0.30 to T = 0.84 and the system size L = 16. 



is given by 

Q3 ^ 2(Q3-'-'-=)3+3g3.^-r= [(g. 



-6Q 



3z — r" 



where is defined in Eq. ( fTSl l. the third order invariant 



+2Q^2^Q^^Q2'^) . (33) 
In a three-dimensional uniaxial nematic state, such as that 
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FIG. 8: (Color online) Spin correlations spanning the nematic state 
and high temperature paramagnet plotted as a function of the rescaled 
distance r\fT on (a) log-log scale and (b) log-linear scales. The 
data are identical to those plotted in Figs. |7Ib) and (c); grey lines 
are guides for the eye, showing ~ l/r** and ~ exp(— r), respec- 
tively. Temperatures should be compared with the nematic ordering 
temperature Tq ~ 0.5. 



realized here, all quadrupoles moments Q" are proportional 
to a simple scalar Q. The presence of a cubic term in the 
free energy Eq. ( l32b therefore implies that the phase transition 
from nematic phase to paramagnet as a function of tempera- 
ture must be first order — as observed in the MC results. 

In principle, thermal fluctuations might select a single Neel 
state from the ice manifold, in which case Tq would mark the 
onset of dipolar as well as quadrupolar order However this is 
not the case. Spin-spin correlations defined in Eq. (|5]i, 

remain short ranged [Fig. |3a)]. At the distances accessible 
to simulation, they rapidly cross over from the power-law 
decay characteristic of the ice manifold at low temperatures 
[Fig. |3b)] to the exponential decay expected for a paramag- 
net [Fig. Etc)]. 

The reason that the usual order-from-disorder mechanism 
is ineffective in selecting dipolar order is the massive degen- 
eracy of the ice manifold — the entropy gain of fluctuations 
about the favored state (relative to the average) would have to 



compensate for the loss of an extensive entropy of 
lnl7o/A^ ^ 0.5 In 1.5 w 0.20 

per spin. We return to this point below. 

At finite temperature, the algebraic decay of spin corre- 
lations D(t) ~ l/r"^ in the Coulomb phase is expected 
to crossover to exponential decay D(r) ^ cxp(— r/^^) for 
r > ^c, where the characteristic length scale diverges 
for r — > 0. For Heisenberg spins in three dimensions, 

~ I/a/ T^'i'"'^'^^ . This is the only length scale in the sim- 
plest Coulomb theory, and it is therefore interesting to plot 
the spin correlations for a rescaled distance r\/T . This 

is done in Fig. |8] At low temperatures T < Tq, the data 
appear to collapse onto a single power-law behavior in this 
range, while they collapse onto an exponential behavior above 
Tq: There is a rapid change between these behaviors, associ- 
ated with the discontinuous transition at T = Tg. The results 
suggest that the Coulomb-phase theory applies to the present 
bilinear-biquadratic model, and in addition, that the character- 
istic length suddenly changes from several lattice spacings 
in the nematic phase for T < Tq to one comparable to the 
lattice spacing in the paramagnetic phase for T > Tq. 

Once again, these results have a simple interpretation in 
terms of local, preformed order In Fig. |6tc) we plot the ex- 
pectation value of the order parameter for the simplest kind of 
four-sublattice order, X^°'°^^ [Eq. ^]. This clearly scales to 
zero with system size. However there is a sharp feature in the 
susceptibility associated with A'^"^^' [Eq. ( fTOb l at Tq, shown 
in Fig. 13 d), where tetrahedra with local E symmetry collec- 
tively choose collinear configurations. Indeed, as T 0, 
A'^*^^' takes on its maximum allowed value of A'^^^' — > 16/3 
[Fig lStc)], as required for loops of perfectly collinear spins. 

B. Nematic phase with local E symmetry 

In applied magnetic field, the "up" and "down" spins of 
the collinear nematic phase immediately "flop" into the plane 
parallel to h, and transform into the 2:2 canted coplanar con- 
figurations shown in Fig. [3] Such canting is entirely com- 
patible with the ice manifold, as is illustrated in Fig. Uc) — 
entire loops of spins cant simultaneously, to give a state with 
smoothly evolving magnetization, but no Neel order. 

The correlation function P(r) retains a finite (reduced) 
value in this new canted manifold of states. However spin 
rotational symmetry is now explicitly broken by the magnetic 
field, so this does not of itself imply nematic order Nematic 
order is none the less present, in the selection of a common 
plane within which the spins cant. This is equivalent to the se- 
lection of a direction (but not an orientation) in the xy plane, 
and long range order can now be observed in the transverse 
moment Q-'--^ defined by Eq. (l22l i. as discussed in TableU 

Since this director breaks the residual 0{2) symmetry, 
the resulting nematic state must possess a branch of gapless 
(Goldstone) modes associated with rotations of the plane of 
canting about the z axis. It is worth noting that a canted Neel 
state with E-type symmetry would break rotational symmetry 
in the same wayi^. However for T-0~^'^ 0, simulations show 
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FIG. 9: (Color online) (a) Nematic order parameter Q^'^ = IQ^'^j 
[Eq. (I22H. signaling the nematic long-range order in finite magnetic 
field; (b) Heat capacity [Eq. ([SB]- Data are at /i = 1.2 and h = 2.4 
for 6 = 0.1, and for system size ranging from L = 4 to L = 16. 

that spin-spin correlations retain their power-law character at 
low temperatures, implying the absence of long-range Neel 
order. 

The transition from local-E symmetry nematic state to 
collinear state as /i ^ is completely smooth, and the finite 
T properties of nematic state at finite h are qualitatively iden- 
tical to those shown in Figs. |6] and |7] with the obvious caveats 
that P(y) < 1 for T ^ 0, and the collinear order parameter 
Q must be replaced by Q-^'^ ~ |Q-'-'^|. Once again the on- 
set of nematic order at Tq is associated with a sharp peak in 
heat capacity, a rise in the local fluctuations with E symmetry, 
Ag*^^', and the absence of long range order of the form A|'°'^^'. 

A suitable free energy to describe this nematic state is 

= a2|Q^''P + C22|Q^''|' + e222|Q-"^'r (34) 

which permits both first and second order phase transitions 
into a paramagnetic phase as a function of temperature, de- 
pending on the sign of C22. However MC simulations sug- 
gest that the transition remains first order. Figure |9] shows 
the temperature dependences of the nematic order parameter 
Q-L^2 ^ |Q_L,2| ^] ^jjjj jjjg jjg^j capacity [Eq. (EBl for 
h ~ 1.2 and h = 2.4. For both cases, the order parameter 
exhibits a sharp onset and the heat capacity shows a jump, in- 
dicating that the E-symmetry nematic transition is of the first 
order, as for /i = in Fig. |6] As noted by comparing the 
results for h = 1.2 and h = 2.4, the discontinuity becomes 
clearer as h increases. 
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FIG. 10: (Color online) Dependence of the magnetization m 
[Eq. 0] on magnetic field h for b = 0.1 and J3 — 0, showing the 
existence of the magnetization plateau in the absence of long-range 
Neel order. Symbols show the result of Monte Carlo simulations for 
temperatures ranging from T — 0.04 to T = 0.24 and the system 
size L — 8. The dashed line is the result obtained by minimizing the 
energy for T — 0. 

In principle the E-symmetry nematic state could interpo- 
late to saturation, simply by canting all spins until they are 
aligned with the magnetic field. However this is not energeti- 
cally favorable at the level of a single tetrahedron (Fig.|3]l, and 
for a magnetic field ft, « 3, the system undergoes a first order 
transition into a state with magnetization m = 1/2, seen as 
the plateau in Fig. [TO] This state is discussed in detail in the 
section below. 



C. Plateau liquid with local T2 symmetry 

The half-magnetization plateaux observed in Cr spinels are 
associated with collinear states with three-up and one-down 
spin per tetrahedron. There are in fact an extensive number 

Qo ~ 1.7^/^ ~ 1.14^ (35) 

of such uuud states — a manifold isomorphic to hard-core 
dimer coverings of the diamond lattice formed by joining the 
centers of tetrahedra^- . (Dimers on bonds of the diamond lat- 
tice correspond to the down spins in uuud states on the py- 
rochlore lattice.) We therefore refer to it as the "dimer" man- 
ifold below. Collinear uuud states with and without simple 
Neel order are illustrated in Fig.fTTI 

It is possible to construct a field theory for the dimer man- 
ifold at T = by exact analogy with the treatment of the 
ice manifold above. The condition that every tetrahedron has 
exactly three-up and exactly one-down spins translates into 
a zero divergence condition for an electric (magnetic) field, 
modified to include a source termii. 

Once again, thermal fluctuations are ineffective in restor- 
ing long-range Neel order. Reduced spin-spin correlations 
Eq. (|5]l at finite distance exhibit a crossover between a dipo- 
lar form [cf. Eq. ( l29] ll for low temperatures, and exponential 
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FIG. 11: (Color online) Half-magnetization plateau states {uuud 
states) on a pyrochlore lattice with exactly three-up and one-down 
spins per tetrahedron, (a) Schematic picture of an uuud state with 
no long range order, associated with the "dimer" manifold; (b) uuud 
state with long-range four-sublattice order with T2 symmetry in- 
duced by FM J3, as considered in Refs. j^] and fisll : (c) 16- 
sublattice order. See also Fig.ll5l 



decay for high temperatures — see Fig. [T2] While spins are 
perfectly collinear at low temperatures [Fig.[T4l'a)l. the z axis 
is now singled out by magnetic field, and Q ^ contributes 
toP(r). 

This means that there is no symmetry breaking associated 
with the smooth rise in collinearity for T* k, b, which should 
be regarded as a crossover rather than a phase transition. Sin- 




FIG. 12: (Color online) Absence of long-range magnetic order in the 
plateau liquid state for b — 0.6, J3 = 0, and h — 4. (a) The reduced 
spin congelation function D{r), defined by Eq. measured along 
the (110) chains, as in Fig.|7ta). (b) Characteristic 1/r^ power-law 
decay of spin correlations at low temperatures T < T* ~ 0.6, asso- 
ciated with the dimer manifold of uuud states, plotted on a log-log 
scale, (c) Exponential decay of spin correlations at high temperatures 
T > T* ^ 0.6, plotted on a log-linear scale. In (b) and (c), grey 
lines show guides for the eye, ~ and ~ exp(— r), respectively. 
In all figures the data are for temperatures ranging from T = 0.30 to 
T = 0.84 and the system size L = 16. 



gular features are similarly absent from the heat capacity, 
shown in Fig. [T4l" b). This smooth change is also seen in the 
rescaled plot of the spin correlations shown in Fig. [13] The 
crossover from the low-T power-law behavior to the high-T 
exponential decay is much more smooth compared to the case 
for the nematic transition at /i = in Fig. |8] This suggests a 
smooth growth of the characteristic length scale ~ 
at T T* « &. We therefore conclude that the magnetiza- 
tion plateau is a true spin-liquid state, continuously connected 
with the high-T paramagnet. We refer to this as the plateau 
liquid below. 
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FIG. 13: (Color online) Spin correlations spanning the plateau liquid 
and high temperature paramagnetic phases, plotted as a function of 
the rescaled distance r^/T on (a) log-log scale and (b) log-linear 
scales. The data are identical to those plotted in Figs.[T2lb) and (c); 
grey lines are guides for the eye, showing ~ l/r'"* and ~ cxp(— r), 
respectively. Temperatures should be compared with the crossover 
scale T* w 0.6. 



It is interesting to note that, despite the absence of any kind 
of long range order, the defining property of the plateau liquid 
— its magnetization (Fig. [TOt — is almost indistinguishable 
from those of the four-sublattice ordered state'^. Long-range 
four-sublattice order is explicitly absent — the plateau 

liquid possess the full Ai symmetry of the paramagnet. None 
the less there is a marked rise in local T2 order of individual 
tetrahedra Aj"' at T* « 6 accompanied by a broad peak in its 
susceptibility, as shown in Figs. [T4|c) and (d). This curious 
spin-liquid state clearly deserves further study. 

To this end, we have performed low-T expansions of the 
free energy of many different ordered and disordered uuud 
states. These are controlled expansions about the ground state 
in powers of T for a spin of length S — \, where we write the 
energy 
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FIG. 14: (Color online) Temperature dependence of (a) collinear- 
ity [Eq. l fT9t . cf. Fig.O, (b) heat capacity [Eq. ||37]|1. (c) the related 
measure of local correlation Aj"' defined by Eq. l llOb . and the global 
order parameter defined by Eq. ((Sj, and (d) the associated lo- 

cal susceptibility. Simulations were performed for /i = 4, b = 0.6, 
in clusters with 1/ = 4 to L = 16. 



in terms of the fluctuations 



(5S = {8Sl,8Sl, . . . 6S'^, SSf, SSl . . . SS^) 



(37) 



^ about a given uuud configuration. The leading fluctuation 

H = Eq ^ — SSiA4ij6Sj + . . . , (36) contribution to the free energy can then be calculated in terms 

^ j J of the trace over eigenvalues of the 2N x 2N matrix Ai in the 
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TABLE II: Fluctuation entropy per site calculated for randomly gen- 
erated uuud states in an A*' = 1024 cluster. Here sf"'", s^""*, and (sf) 
are the lowest, highest, and mean value of the entropy, respectively, 
and Asf = s^"'^ — (sf) measures the deviation of the highest value 
of entropy from the mean. Statistical errors on all numbers are less 
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(38) 



where Eo is the ground state energy, its degeneracy, and 
(...)no the average over all degenerate ground states. 

For a generic ordered phase, ^Iq is finite, and del A4 takes 
on the same value for all (symmetry related) ground states. In 
this case Inflo/N — )• for N ^ oo. However for the dimer 
manifold, Oq « 1.14^, which means that the ground state has 
a finite entropy per site 



So 
N 



In 1.14 « 0.13 



(39) 



In this case, different ground states are not related by simple 
lattice symmetries and the fluctuation entropy per site 



Sf 



In dot M 
2N 



(40) 



takes on a range of values. 

We have studied the distribution of values of Sf within 
the dimer manifold for a range of values of b, by numeri- 
cally calculating dot Ai for 10000 randomly generated uuud 
states in a cluster of = 1024 sites (L = 4), using a 
Monte Carlo algorithm based on loop updates of spins. We 
found that the highest value of Sf is achieved by an eight- 
fold degenerate, 16-sublattice "R-state"— , in which the four 
A-sublattice tetrahedra within the 16-site cubic unit cell of 
the pyrochlore lattice take on all four possible uuud config- 
urations [Fig. [TTlc)]. This state has overall cubic symmetry, 
and is actually observed in the plateau phase of HgCd204^. 
The lowest value of Sf is achieved by the four-sublattice or- 
der shown in Fig. [TTT b). The calculated values of the maxi- 
mum and minimum values s^^'^ and sf"'" are listed in Table HIl 
together with the mean value (sf ) and the difference between 
sf"^^ and the mean (sf), Asf . 

From these results it is immediately clear why thermal fluc- 
tuations alone fail to select a unique ground state for any value 



FIG. 15: (Color online) (a) Probability distribution of the flippable 
hexagons within the dimer manifold, (b) The fluctuation entropy per 
site Sf [Eq. J40b 1 as a function of the number of "flippable" hexagons. 
The lower bound Sf = —4.19874 is set by the four-sublattice state 
shown in Fig. ll If b). which has no flippable hexagons. The upper 
bound Sf = —4.19806 is set by the 16-sublattice state with the 
maximum number of flippable hexagons [see Fig. ll lf c)1. The blue 
dots represent a sample of 10000 random configurations, (c) Prob- 
ability distribution of the fluctuation entropy per site Sf within the 
dimer manifold. All results are for a cluster of A" = 1024 sites with 
6 = 0.6. 



of b considered in this paper The fluctuation entropy per site 
gained by choosing the cubic 16-sublattice state is miserly, 
for example, Asf = 0.00183 for b = 0.1 and Asf = 0.00051 
for b = 0.6. These numbers must be compared with the ex- 
tensive entropy Sq/N w 0.13 of the liquid phase, all of which 
is lost if the system orders. So for the values of b considered 
here, thermal fluctuations cannot drive the system to order. 

However it is amusing to note that the entropy gain Asf in- 
creases as b decreases, scaling approximately as In b, as shown 
in Tablellll This raises the intriguing possibility that b acts as a 
singular perturbation, and that for sufficiently small b, fluctu- 
ations might overcome the extensive entropy Sq/N « 0.13 of 
the dimer manifold, driving the system order — even though 
it is disordered for 5 = 0. Such an order-from disorder ef- 
fect would presumably favor the cubic 16-sublattice R-state, 
which is also believed to be selected by quantum fluctuations 
at T = q12j13j37 However in the present model, it would occur 
only for vanishingly low temperatures, and would therefore 
be extremely difficult to access in simulation. This question 
remains for future study. 

The result above explains why the system does not order 
at finite temperature, but not why the fluctuation entropy fa- 
vors the 16-sublattice state? We can answer this question 
by looking at the distribution of the fluctuation entropies Sf 
within the dimer manifold of uuud states. Figure [TSl shows 
the distribution for 6 = 0.6. The uuud states can be broken 
up into classes of states with a different net flux of an effec- 
tive magnetic (or, equivalently, electric) fiel d'^-'-^i^^'-^'i^^ . This 
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flux is conserved by all cyclic exchanges of spins on loops of 
alternating u and d spins, motivating a loop expansion of the 
fluctuation contribution to the free energy^^. The leading 
term in such an expansion counts the number A'fiip of six-site 
hexagonal rings in a "flippable" u-d-u-d-u-d configuration. 

In Fig. [TSlb) we plot the fluctuation entropy Sf as a func- 
tion of iVfijp. The highest (lowest) values are achieved for the 
16-sublattice (four-sublattice) ordered states with the most 
(least) flippable hexagons. The fact that randomly generated 
uuud states lie extremely close to the line connecting these 
two states suggests that loops of more than six sites contribute 
little to the fluctuation entropy. 

The remaining question is why the overall difference in 
fluctuation entropy between different uuud states is so small? 
We can express the fluctuation entropy in terms of the 2N 
eigenvalues {e„} of M. as 



Sf = 



2N 



(41) 



The eigenvalue spectrum {e„} associated with the simplest 
q = four-sublattice uuud state can easily be calculated an- 
alytically; working with a four-site unit cell, there are four 
bands, one of which is nondispersive. The associated density 
of states (DOS) for b = 0.6 is shown in Fig. [T6l a). where 
the flat band appears as a sharp peak at e = 166 + /i — 4. 
In Figs. [T6l b)-(d) we compare the integrated DOS from these 
four bands with numerical results for the integrated DOS of a 
1024— site cluster. 

The integrated DOS, averaged within the dimer manifold of 
uuud states [Fig.[T6l'd)l. is indistinguishable by eye from that 
of the four-sublattice state [Fig. [T6l b)l. The step associated 
with the flat band survives as a set of N /A localized excitations 
at e = 165. And, critically, the gap 



A 



i + 4- V(85 + 4)(8& + 4-/i) + /i2 



(42) 



to the lowest lying excitation is set by a nodeless eigenvec- 
tor, whose components SSf depend only on whether the spin 
Si points up or down. All uuud states can be made formally 
equivalent to four-sublattice order by renumbering the sites 
in each individual tetrahedron, and the energy of this nodeless 
excitation is also unchanged by this renumbering of sites. It is 
therefore completely insensitive to whether or not the system 
is ordered. From the results it is clear why the thermody- 
namic properties of the plateau liquid state, and in particular 
the entropy associated with fluctuations about it, are so close 
to those of the ordered plateau state. 

From these results, it is also possible to understand why the 
numerically determined entropy gain Asf increases as 6 
for = 4 (cf Table HHi. This singular behavior can be traced 
back to a band of excitations above the spin- wave gap A w 4b, 
with bandwidth Ae ^ b, which collapses to become a strict set 
of zero modes for 6 — > 0. Since zero modes are excluded from 
the sum which determines Asf , while the collapsing band con- 
tributes as ~ In 6, b acts as a singular perturbation, and in- 
finitesimal b may drive the system to order. This is despite the 
fact that it is disordered for = 0, and for the relatively large 
b used in our simulations. 
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FIG. 16: (Color online) (a) Density of states (DOS) for eigenvalues 
of A4 for four-sublattice uuud state in thermodynamic limit, show- 
ing finite gap A [Eq. ( I42H and flat band at e = 16b. Integrated DOS 
(IDOS) (b) of four-sublattice uuud state (points), (c) of a typical 
disordered uuud state (points), and (d) averaged within disordered 
uuud states (dashed line). In all cases Js — 0, h = 4, and b = 0.6, 
and a cluster size is A*" = 1024. The integrated DOS corresponding 
to (a) is shown in shading on (b)-(d) for comparison. 



D. vector-multipole phase with local T2 symmetry 

At the upper critical field of the magnetization plateau, the 
collinear spins of the uuud configurations cant away from the 
z axis. This instability occurs at the level of a single tetrahe- 
dron (Fig.O, where it is continuous. On a lattice, it is associ- 
ated with the closing of the gap A [Eq. ( l42b l in the excitation 
spectrum of the plateau liquid. Because of the special struc- 
ture of this excitation, discussed above, the gap closes at the 
same value of /i^ = 4 + 86 for all uuud states, and the transi- 
tion is once again continuous — at least for T = 0. However, 
since the spin configurations in question are simply 3:1 canted 
versions of the uuud states, with local T2 symmetry, all of the 
entropic arguments presented above for the plateau liquid still 
hold. Thermal fluctuations alone cannot restore (canted) Neel 
order, and spin-spin correlations exhibit a power-law decay 
of l/r^ forT 0. 

The resulting state does however exhibit long range order 
in both the rank-two tensor order parameters Q-'-^^ and Q-'-'^ 
[Eqs. ( l22l i and ( l23T l. and TableU). The 3:1 canting of the uuud 
spins selects a direction in the xy plane, and the primary order 
parameter is therefore the lower-symmetry irrep, Q-'--^. The 
finite value of the nematic order parameter Q-'-'^ reflects the 
fact that this canting is coplanar Since Q-^'^ transforms like 
a vector under rotations about the z axis, we classify this state 
as a vector-multipole phase with local T2 symmetry. 

Within the framework of a Ginzburg-Landau theory, the 
contribution to the free energy from this pair of order param- 
eters is 
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FIG. 17: (Color online) Magnetic field dependence of (a) primary 
order parameter Q^'^ = [Eq- i23b 1 and (b) secondary order 

parameter Q^'^ = |Q^'^| [Eq. j22M in the vector-multipole phase 
for 6 = 0.1. The continuous transition from the plateau liquid state 
into the vector-multipole phase at lower magnetic fields, and the di- 
rect transition from the vector-multipole phase into the (saturated) 
paramagnet at higher fields are clearly visible. Dashed lines show 
behavior at T = in a single-tetrahedron theory. Points show results 
of MC simulations for the system size N = 8 from T — 0.02 to 
r = 0.08. 



and Q^'^ 



and (Ell), 



-L,2 



where, following Eqs. 

Q2 ^ g^^i Qi 

Eq. ( l43T l should be contrasted with the form of free en- 
ergy in the absence of magnetic field, Eq. ( |32] |. The cubic 
invariant Q"^ survives as an interaction 612 between Q-'--^ and 
Q-'-'-^, which transforms like e2j0g-i0g-i0 ^ ^ under rota- 
tions about the z axis. This means that components of Q-'-'^ 
couple linearly to a quadratic combination of the components 
of Q-'-'^. Because of this, a finite value of the (lower symme- 
try) order parameter Q^'^, immediately induces a finite value 
of the (higher symmetry) order parameter Q-'-'^. 

In principle Eq. ( |43] ) permits both first and second order 
phase transitions into the vector-multipole phase from disor- 
dered (paramagnetic) or pure nematic phases, depending on 
the sign of the coefficients cn, C12, and 022- The full solution 
for Q-'-^^ and Q-'-'^ is further complicated by the fact that these 
order parameters also couple to octupolar spin moments (see 
Appendix iBl for details). However the relationship between 
Q-*-'^ and Q-'-'^ is clear at the level of a single-tetrahedron 
theory (cf. Ref. lH). 

Within the theory for a single, embedded tetrahedron — 
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FIG. 18: (Color online) Temperature dependence of (a) the primary 
order parameter Q^'^ = |Q^'^| defined by Eq. l l23b . showing the 
onset of the vector-multipole order at T = Tv ~ 0.042, (b) heat 
capacity Eq. Pit , (c) the related measure of local correlation Aj"' 
defined by Eq. dlOt , and the global order parameter defined by 
Eq. {8](, and (d) the associated local susceptibility Eq. i ll It . Simula- 
tions were performed for /i = 6, 6 = 0.1, in clusters with 1/ = 4 to 
L = 16. 



which is exact for T — — the primary order parameter 

Q-^'^ = IQ''"'^! grows as 



: \/h — he 



(44) 



^2(3 - 2b) 

while the secondary order parameter Q-^'^ — IQ^'^I grows 
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more slowly as 

Q^-'ih>h,)^^^^^^ih-K). (45) 

The results of this theory for the T2 vector-multipole phase 
are shown by the dashed lines in Fig. [17] For the value of 
b used in the present study, the zero temperature transition 
from T2 vector-multipole phase to paramagnet at high field 
is strongly first order, even at the level of a single-tetrahedron 
theory — cf. Fig.|3] — and remains so throughout. 

The nature of the finite temperature transition from the 
vector-multipole phase into the paramagnet is harder to de- 
termine. However, as shown in Fig. [17] it appears to be first 
order for all h> h*, where (T*, /i*) « (0.1, 5.2) marks the 
point at which the crossover line T* joins the boundary of the 
vector-multipole phase, Ty, as shown in Fig. |2tb). On the 
basis of our results, we consider that there is a tricritical point 
at (T* , h* ) where the nature of the phase transition into the 
T2 vector-multipole phase changes from continuous to first 
order. 

In Fig. [18] we present MC simulation results for the finite- 
temperature transition into the vector-multipole phase for 
h = 6.0 > h*. The primary order parameter Q-'-'^ be- 
comes nonzero with a sharp jump at a transition temperature 
Tv ~ 0.042 [Fig. [TSa)]. Both heat capacity and local T2 
susceptibility show a jump at the transition Ty [Figs. [TSl b) 
and (d)], but no sign of long range order in the bond-order 
parameter given by Eq. ((TJ [Fig.[T8jc)]. 

The finite temperature transition from the plateau liquid to 
T2 vector-multipole phase for h < h* deserves special atten- 
tion, since the plateau liquid exhibits algebraic decay of cor- 
relations for intermediate distances. Monte Carlo simulations 
suggest that the transition has a continuous character with (ap- 
proximately) mean-field exponents. We return to this below 
in Sec HV^ 



E. Global structure of the h-T phase diagram 

Our results for the h-T phase diagram of the antiferromag- 
netic nearest-neighbor Heisenberg model with additional bi- 
quadratic interactions b [Eq. ([TJ] are summarized in Fig.|2jb). 
There are two ordered phases, a nematic phase with local E 
symmetry and a vector-multipole phase with local T2 sym- 
metry, both of which break spin rotational symmetry about 
the direction of the magnetic field. These are separated by a 
plateau-liquid state with all the symmetries of a paramagnet 
in magnetic field. 

This phase diagram bears a very strong resemblance to that 
of the corresponding model with weak FM third-neighbor in- 
teraction J3 ~ —0.05, which enforces four-sublattice order, 
as shown in Fig. |2a) (cf. Ref. ITsll ). So far as the topology 
of the phase diagram is concerned the only change is the re- 
placement of a line of first order phase transitions terminating 
the four-sublattice plateau state (which breaks lattice symme- 
tries), by a crossover in the case of the plateau liquid (which 
does not). 



Throughout this paper, we have argued that preformed lo- 
cal order at the level of a single tetrahedron exists in all of 
these phases. Moreover, in the case of the half-magnetization 
plateau, we have seen in Sec. 1111 Ci that conventional magnetic 
order has very little impact on the excitation spectrum, and 
therefore on the thermodynamic properties of the system. 

Viewed in this way, the correspondence between the two 
h-T phase diagrams is not at all surprising — the role of sec- 
ondary interactions like J3 is merely to select between an in- 
finite set of different ordered ground states. Precisely how 
enforcement of long range order works at finite temperature is 
a complex and very interesting question, to which we provide 
only a partial answer below. 



IV. THERMAL TRANSITIONS BETWEEN DIFFERENT 
ORDERED AND DISORDERED STATES 

A. General context 

None of the phases described above possess conventional 
magnetic order of the form (Si) 7^ 0. However they all con- 
tain the seeds of such order in the form of well formed local 
orders A^"' and A!^^"'. Long range order can easily be re- 
stored by adding additional terms to the Hamiltonian Eq. ([TJ. 
The simplest possible choice is a FM third-neighbor interac- 
tion J3 < in Eq. (O, leading to four-sublattice order of the 
form considered in Refs. 10] and ITsll . In what follows we 
study how FM | J3I ^ b precipitates an ordered uuud state 
from the plateau liquid for h ~ A, and contrast this with the 
way in which Neel order emerges from the nematic phase with 
local E symmetry for /i = 0. We also discuss the continuous 
transition from plateau liquid to vector-multipole phase for 
h < h*, T < T*. 

We study these phase transitions as a function of tempera- 
ture T which also gives us access to the high temperature para- 
magnetic phase. This is interesting because, for intermediate 
distances r < ^ VT the nematic phases exhibit the alge- 
braic decay of spin correlations characteristic of a Coulomb 
phase, rather than the exponential decay of correlations more 
usually associated with a paramagnet. Transitions between a 
disordered phase subject to an ice-rule type constraint and a 
phase with conventional order have been discussed for a long 
time in the context of hydrogen-bonded ferroelectrics^^. 
More recently such questions have arisen again in the con- 
text of experiments of many highly frustrated magnets™, and 
in the past few years there has been a theoretical effort to un- 
derstand how order can emerge from a Coulomb phase in clas- 
sical dimer— and spin models^^i^. 

A strong motivation for this work has been the possibility 
of observing an unusual continuous phase transitions, includ- 
ing transitions lying outside the Landau-Ginzburg-Wilson 
paradigm^^. Indeed the (classical) dimer model on cubic 
lattice does exhibit a continuous transition from a Coulomb 
phase at high temperatures to a simple crystalline ordered 
phase as a function of temperature'*! This transition has 
unusual scaling properties-"^, and does not naively admit a 
Landau-Ginzburg-Wilson description, since the high temper- 



17 



ature phase cannot be described using an expansion in terms 
of the low-temperature order parameter A recent very de- 
tailed simulation study of a family of three-dimensional dimer 
models with ordered ground states and high-temperature 
Coulomb phases found a rich variety of continuous and dis- 
continuous phase transitions, including double phase tran- 
sitions where monopole excitations condense out of the 
Coulomb phase to give a conventional paramagnet at inter- 
mediate temperatures^a. 

The present understanding of these phenomena is that the 
gauge field associated with Coulomb phase is minimally cou- 
pled to a matter field which condenses in the ordered phase, 
following an Anderson-Higgs mechanism^— . In fact it is 
also possible to study zero temperature (quantum) phase tran- 
sitions from Coulomb to ordered phases in three-dimensional 
quantum dimer modelsSS. These can in principle be continu- 
ous, occurring through the condensation of monopole excita- 
tions in the Coulomb phase--, but numerical simulations sug- 
gest that the transition is first order^*^. 

Less is known about transitions in spin models, but one 
interesting scenario exists for a continuous transition in an 
extended Heisenberg model on a pyrochlore lattice from a 
Coulomb phase to a four-sublattice ordered stated. This tran- 
sition is found to be in the same universality class as a uniax- 
ial ferroelectric with dipolar interactions, for which the upper 
critical dimension is three^*^. This makes possible to continu- 
ous transitions with mean-field exponents (up to log correc- 
tions) — a scenario which closely resembles the transition 
from plateau-liquid into vector-quadrupole phase discussed 
below. Generically, however, transitions from Coulomb liq- 
uids into ordered states seem to be first order~, a fact which 
may be explained by interactions between fluctuations of as- 
sociated gauge field^. 

We conclude by noting that the complex forms of order 
which can occur in Heisenberg models on the pyrochlore lat- 
tice as a result of the interplay between farther-neighbor in- 
teractions and thermal fluctuations are also a topic of cur- 
rent interest^^. In finite magnetic field, these lead to a half- 
magnetization plateau which can be tuned at will between dif- 
ferent forms of order—. A similar fluctuation driven plateau, 
but with a uniquely defined form of order, is also expected to 
occur for the edge sharing tetrahedra of the FCC lattice^. 

We now return to the model in question. 



B. Transition from plateau-liquid to ordered uuud state 

For h ~ A, T < h, Eq. (dJ exhibits the plateau-liquid 
state described in Sec. 1111 CI [cf. Fig. |3b)]. Inclusion of a 
FM third-neighbor interaction J3 [Eq. (|2|i] causes it to or- 
der at low temperatures. We consider first the conventional 
limit where both | J3I and h are "large", choosing parameters 
J3 = —0.06 and h ~ 0.6. In this case there is strongly first or- 
der transition from paramagnet to four-sublattice plateau state 
for T/v ~ 0.70. This can be seen very clearly in simulation re- 
sults for the heat capacity and the order parameter Aj°'^^', and 

its susceptibility Xj°'^^', presented in Fig. [191 If we now de- 
crease I J3I, the transition temperature Tm must also decrease. 
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FIG. 19: (Color online) Temperature dependence of (a) heat capacity 
defined by Eqs. (jSTJ, (b) and (c) the related measure of local corre- 
lation Ay"' and its susceptibility defined by Eqs. l llOt and jilt , and 
(d) and (e) the global order parameter and its susceptibility de- 
fined by Eqs. ^ and l|9} for h = 0.6, h = A and a range of values of 
system size. The results are for J3 = —0.06, showing a single first 
order transition into the ordered phase at Tjv = 0.70(1) (indicated 
by the vertical dashed line). 



and for sufficiently small | J3 1 it will become smaller than the 
crossover temperature T* w 6 associated with the plateau liq- 
uid. 

In this case, there are anomalies in thermodynamic quanti- 
ties at two distinct temperatures as demonstrated in Fig. |20] 
There is a broad maximum in Xj""' at T* « 0.6 [Fig.|Mc)], 
signaling the onset of the plateau liquid state, accompanied by 
a broad peak in the heat capacity at a slightly lower temper- 
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FIG. 20: (Color online) The same plots as Fig. [T9] but for 
Jg = —0.02, showing both a crossover in to the plateau liquid for 
T* ~ 0.6 (bold grey line) and a transition into the ordered phase at 
Tjv = 0.46(1) (dashed line). 




FIG. 2 1 : (Color online) Temperature dependence of the order param- 
eter for four-sublattice uuud order [as defined by Eq. Isjl], for 
a range of values of J3 < 0. The transition temperature becomes 
smaller as j J3 1 — >■ 0. At the same time the transition becomes more 
strongly first-order. All results are for h = A and h = 0.6, in a cluster 
with L = 8. The lines are guides for the eye. 
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FIG. 22: (Color online) Phase diagram for the classical Heisenberg 
antiferromagnet on a pyrochlore lattice in applied magnetic field 
ft = 4, with additional biquadratic interactions b — 0.6 [Eq. ([TJ]. 
The transition temperature Tjv associated with the gapped, ordered, 
half-magnetization plateau state vanishes as the strength of ferro- 
magnetic third-neighbor interactions J3 — >■ 0, as determined by 
Monte Carlo simulation. A state exhibiting a half-magnetization 
plateau but no long-range magnetic order exists above Tn up to a 
crossover temperature T*. Estimates of the crossover temperature 
are taken from peaks in the local susceptibility and heat capacity. 



ature [Fig.lMa)]. And, at T = T^r « 0.46 < T*, there is a 
small jump in Cy, accompanied by a clear singularity in global 
order parameter susceptibility Xjf'^^ [Fig.l20l'e)l. While there 
is no true phase transition at T*, it is clear that the bulk of the 
entropy of the paramagnet is lost in the smooth crossover into 
the plateau liquid, and not in the first order transition into the 
ordered phase. 

So what happens for J3 — > 0? Unfortunately this question 
is hard to answer by Monte Carlo simulation, as the massive 
degeneracy of the uuud states translates into many competing 
local minima in the free energy. However T/y is strictly zero 
for J3 — 0, and there are two obvious scenarios for how this 



can be achieved. 

The first is that the first-order transition into the ordered 
phase becomes weaker as Tjv 0, terminating in a critical 
end point for J3 = 0, Tjv = 0. This end point would in fact 
be multicritical, since many different ordered uuud states can 
be formed out of the dimer manifold for different choice of 
long range interactions. Within this scenario, the power-law 
correlations between spins in the plateau liquid for T — > 
could be viewed as evidence of critical fluctuations. The sec- 
ond scenario is that the first order transition into the ordered 
phase persists down to T/v = 0. Since an infinite number of 
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different ordered phases branch out from the point J3 = 0, 
Tjv — 0, it can probably best be termed "multifurcative" . 

First-order phase transitions between different ordered 
phases with an infinite degeneracy at the transition occur in 
a number of models. Such phase transitions are first order, in 
the sense that neither ordered parameter collapses approach- 
ing the critical point. However they also exhibit one of the 
characteristic features of a second order transition, namely a 
soft excitation or set of soft excitations connecting the differ- 
ent ordered ground states. 

As far as we can tell from our present results, it seems most 
likely that the classical pyrochlore AF with biquadratic in- 
teractions exists at a multifurcative point in parameter space, 
with an infinite ground-state degeneracy, not at a critical end 
point. As shown in Fig.|2Tl the low temperature value of the 
order parameter is broadly independent of Tn. This implies 
that the phase transition in fact becomes more strongly first 
order at T^v ^ 0, and appears to rule out a (multi)critical end 
point. Our collected simulation results for J3 are sum- 
marized in the form of the phase diagram in Fig. |22] 

It is amusing to note that this phase diagram bears a su- 
perficial resemblance to the phenomenology of a second or- 
der (quantum) critical point — a transition temperature which 
collapses to a special point with algebraic decay of correlation 
functions, which in turn controls a broad region of the phase 
diagram up to a characteristic crossover temperature T* . All 
of this despite the fact that the only phase transition present 
is first order, which means that the length scale associated 
with fluctuations remains finite. Some of the generic features 
seen in our model — power law decay of correlations over 
a large, but finite, length scale — have been previously dis- 
cussed in the context of other models with strong local con- 
straints, where they were dubbed "high temperature critical- 
ity"i2. 

The transition from a critical "Coulombic" phase described 
by a ?7(1) gauge theory into a simple ordered state as a func- 
tion of temperature can be studied much more cleanly in the 
(classical) dimer model on cubic lattice, where the constraint 
enforcing the dimer manifold is infinite. In this case, the phase 
transition is continuous, and exhibits interesting and unusual 
scaling propertie a^^'^^ . We have made a preliminary study of 
the "stiffness" K associated with fluctuations in a U{1) gauge 
theory for temperatures spanning the paramagnet and plateau 
liquid phases in our model (see Fig. l22T i. but find no clear ev- 
idence of a phase transition. However the way in which the 
dimer and loop manifolds break down at finite temperature 
in a model with a finite constraint is an interesting problem, 
and one which deserves further study. We note in passing that 
interesting, related, problems arising the context of quantum 
loop models^. 



C. Transitions from paramagnet to E-symmetry nematic 
pliase and Neel ordered state 

It is interesting to contrast the finite temperature phase tran- 
sitions associated with the plateau states for h w 4, with 
the transitions into E-symmetry Neel and nematic ordered 



states for h = 0. Once again, for "large" IJ3I there is a 
strongly first order transition from the paramagnet into the 
Neel phase at a unique temperature Tjy. Meanwhile, for 
"small" I J3 1 ^ b, there is double transition, first from para- 
magnet to E-symmetry nematic phase at Tq ^ b, and then 
into the four-sublattice Neel order at a much lower tempera- 
ture Tjv, as demonstrated in Fig. [23] Within the limits of our 
simulation, both of these transitions appear to be first order in 
character^^. The results for varying J3 are summarized in the 
phase diagram in Fig.|24] The first order transition from para- 
magnet to nematic phase at Tq at /i = should be compared 
with the crossover from paramagnet to plateau liquid T* ^ b 
observed for /i = 4 in Figl22] 



D. Transition from plateau-liquid to vector-multipole phase 

Perhaps the most interesting of the finite temperature transi- 
tions observed in our model is the one from plateau-liquid to 
vector-multipole phase, already described in Sec. IIIIDI At 
one level this is the most exotic phase transition we study 
— a continuous phase transition from a "Coulombic" state 
with algebraic decay of correlation functions (the plateau liq- 
uid) to a phase with long-range multipolar order (the vector- 
multipole state). But at the same time it has the simplest 
phenomenology of any of the phase transitions in this pa- 
per, with the order parameter exhibiting a simple mean-field 
like behavior Q-^-^{T) - y^Ty - T with Ty ~ 0.09, as 
shown in Fig. l25T a). The secondary order parameter Q^'^ 
grows more slowly as expected [Fig. |25l b)1. The heat ca- 
pacity does not show a noticeable singularity at T = Ty in 
Fig. l25l c). which is also consistent with the mean-field be- 
havior Cy - {T -Tv)" with a = 0. (The broad peak at 
T ~ 0.112 again corresponds to the crossover temperature 
T* for the plateau-liquid state.) 

At a qualitative level, and in the spirit of this paper, it is easy 
to see how a continuous transition can arise between these two 
states. Both are built of tetrahedra with a local T2 character, 
with three "up" and one "down" spin, joined at the corners. 
Both states will exhibit algebraic decay spin correlations at 
low temperatures, as a result of the infinite number of ways 
that these tetrahedra can be assembled to form a pyrochlore 
lattice. The only difference is that three "up" and one "down" 
spins are canted in the vector-multipole phase, giving a finite 
value of Q-^'^ and Q-'-'^ [Figs.|23a) and (b)]. As long as this 
canting can interpolate smoothly to zero in the collinear uuud 
state, the transition will be continuous. And at the level of a 
Ginzburg-Landau theory, nothing prevents this from happen- 
ing — cf. Eq. ( |43] |. 

However in principle it should also be possible to tran- 
scribe each of these phases in terms of the more sophisticated 
"solenoidal field" theory used to describe Neel order in a spin 
model with a high temperature Coulomb phase (cf. Ref. [53[|). 
To the best of our knowledge, nobody has yet attempted to 
extend the gauge-field description of a Heisenberg type spin 
model to treat multipolar order. But it is interesting to note 
that the transition from Coulombic phase to simple Neel or- 
der was found to be continuous, and in a universality class 
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with upper critical dimension three, i.e., one where the critical 
behavior is mean-field like, up to log corrections^^*^. 



V. SUMMARY AND CONCLUSIONS 

We have studied the ordered and disordered phases of the 
classical, bilinear-biquadratic Heisenberg model on the py- 
rochlore lattice at finite temperature and in applied magnetic 
field. We find a rich collection of unconventional states — ne- 
matic and vector-multipole phases with distinct and different 
local symmetries, separated by a half-magnetization plateau 
with spin-liquid character. All of these phases show an under- 
lying "Coulombic" character with algebraic decay of spin cor- 
relation functions over distances ?" < ■Cc ~ 1/ V^- Interest- 
ingly, the transition from plateau-liquid to vector-multipole 
phase is continuous, and appears to be well-described by 
mean field theory. 

While this behavior is undeniably exotic, all of these states 
can be understood — and even anticipated — from a proper 
understanding of the geometry of the pyrochlore lattice, and 
the properties of a single tetrahedron. Strong local fluctua- 
tions of Neel order are present in all of these phases, and 
the zero temperature phase diagram can be understood sim- 
ply from the "self assembly" of these ordered tetrahedra into 
complex states with higher symmetry. 

It is therefore unsurprising that conventional Neel order 
(with four-sublattice structure) is immediately restored by 
the introduction of a ferromagnetic third-neighbor coupling 
J3. However for small | J3I, the unconventional states survive 
above the Neel transition temperature Tn- In particular, the 
spin-liquid plateau survives above Tjv, up to a crossover tem- 
perature T* « b. The transition between liquid and ordered 
plateaux is first order in nature, and remains so for T/v — > 0. 
For small | J3I, the system also exhibits a first-order transition 
between the Neel and nematic phases, in addition to the first- 
order transition from the nematic phase into high-temperature 
paramagnet. 

So far as experiment is concerned, our main finding is that 
the physics of a pyrochlore antiferromagnet in magnetic field 
can be largely determined by the properties of a single tetra- 
hedron. In the simple models which we have considered it is 
possible to tune between states with entirely different point 
group symmetries at will, simply by changing the form of 
(weak) long range interactions present. This is an oversimpli- 
fication, in the sense that magnetostriction in real systems is 
likely single out a particular phonon (or family of phonons) 
with definite symmetry, which will then drive the system 
towards collinearity {b, in our model) and select the low- 
temperature ordering pattern (long range interactions, e.g., J3, 
in our model). However, as long as there is a strong coupling 
to phonons within individual tetrahedra, the form of the mag- 
netization plateau and associated phases may be largely inde- 
pendent of these (system dependent) details. 
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Appendix A: Classification of symmetry breaking at the level of 
a single site 

In order to identify the different possible forms of magnetic 
order which can survive where conventional Neel order breaks 
down, it is helpful to classify the different forms of symmetry 
breaking which exist at the level of a single site. This analysis 
is in the spirit of the detailed classification for the nematics 
in liquid crystals undertaken in Ref. ll62|] . and motivates the 
rank-two and rank-three tensor order parameters introduced 
in Sec. Ill CI and Appendix IB] In order to keep contact with 
quantum spins, which are axial rather than polar vectors, we 
must keep track of time reversal symmetry. 

In Table HU] we show the transformation rules for the spins 
under selected symmetry operations, including time reversal 
0S = — S. In contrast to the usual polar vectors, inver- 
sion leaves the axial vectors invariant - as a consequence, all 
the usual (reflection, rotation, and inversion) symmetry oper- 
ation can be represented by an orthogonal matrix belonging to 
50(3), with determinant equal to H-1. The role of inversion 
in the case of polar vectors is taken over by the time reversal 
operator 6. 

All of the symmetry operations, extended with the time re- 
versal, can be represented by orthogonal matrices with deter- 
minant -1. In the Table HHI we also check if the collinear and 
coplanar states are invariant under those symmetry operations. 
Since we are interested in the symmetry breaking which can 
occur in the absence of broken translational symmetry, we do 
not apply the symmetry elements to the lattice points (i.e., we 
treat all the spins as they were at the origin). We find that the 
invariant operations of the 2:2 state include a C2(z) rotation 
in addition to the symmetry operations of the 3:1 state. 

In Table |IV] we show the symmetry group of each of the 
spin states. In order to facilitate comparison with Ref. ll62ll . 
we also show the symmetry group of the states if the spins 
were polar vectors. We can see that as the time reversal does 
not play a role for the 2:2 collinear state, its symmetry group 
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being the grey-group Dock + ©^co/i- The magnetic (the 3:1 
collinear and both canted states) states have a magnetic point 
group as a symmetry group. 

When studying which symmetry group is broken for the 
magnetic states, we need to note that the external magnetic 
fields lowers the 0(3) symmetry of the space to Coo x 
{E, 1} + QayCoo X {E, I}, where the axis of the Coo is par- 
allel to the magnetic field, and CTu is a reflection to a plane 
that includes the axis of the magnetic field. The symmetry 
of the space with magnetic field is actually identical to the 
symmetry of the 3:1 collinear state. Thus, within a Ginzburg- 
Landau framework we do not expect a continuous (second or- 
der) phase transition between the T = liquid plateau and the 
high temperature disordered phase. The Z2 lowered symmetry 
of the 3:1 state with respect to 2:2 canted state is manifested in 
the Z2 symmetry lowering of the vector to the nematic phase. 



Appendix B: Higher order multipoles 

In this paper, we have classified states according to the low- 
est moment of spins which breaks spin rotational symmetry. 
According to this conventional, "common sense" prescription, 
a state which lacks conventional dipolar (e.g., Neel) order, but 
exhibits a common plane for the canting of spins, is automati- 
cally classified as a nematic or vector-multipole phase. While 
this classification scheme is unambiguous, it is not complete, 
and in some cases may give the wrong answer, so far as the 
primary order parameter is concerned. 

This point was recently discussed at length for the copla- 
nar ground-state manifold of the classical Heisenberg model 
on a kagome lattice, where the primary order parameter was 
convincingly argued to be octupolar, and yiot quadrupolar, in 
nature^i. Incorrect assignment of the primary order parame- 
ter does not affect our ability to detect a bulk ordered phase, 
but can lead to false conclusions about phase transitions. This 
is particularly true of two-dimensional systems at finite tem- 
perature, where the homotopy group associated with the order 



TABLE III: The transformation of spins under different symmetry 
operations. E is the identity element, / is the inversion, O is the 
time reversal operation, 0-^13 is a reflection with a mirror plane q/3, 
and C2(a) is a two-fold rotation around the a axis. In the last two 
columns we indicate if the 2:2 and 3:1 canted states (with magnetic 
moment along the z axis and spins are in the xz plane) are invariant 
with respect to the particular operation. 
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parameter determines the form of topological defect entering 
into Berezinsky-Kosterlitz-Thouless type phase transitions. 

In fact the states which we classify as "nematic" or "vector- 
multipole" in Sec. Hn] also posses higher order multipole mo- 
ments which, under some circumstances, couple to the rank- 
two tensor order parameters used in this paper. We illustrate 
this below for the specific case of the rank-three tensor asso- 
ciated with octupolar order. 

This is odd under time reversal, and has seven components 



^ 1 



given by 

i 

rpV^ 



-3yx 



I 1 i\ 

-V ) _ 



T, 



Sf [(5f )2 



(Bl) 



(B2) 
(B3) 
(B4) 
(B5) 

(B6) 



3{Sff ~ 2(5f )2](B8) 



(B7) 



In the absence of magnetic field, quadrupolar order can cou- 
ple to (fluctuations of) octupolar order through terms of the 
form 



(B9) 



in the free energy, which respect the full 0(3) symmetry of 
the Hamiltonian, and time reversal invarianco^. Therefore, a 
finite octupolar order parameter usually induces a quadrupolar 
one, while the opposite is not always true. When they occur 
together, some care must then be taken to assign the correct 
primary order parameter 

In finite magnetic field we again classify these octupoles 
according to the way in which they transform under rotations 
about direction of magnetic field (the z axis). We obtain a sin- 
gle one-dimensional irrep and three two-dimensional irreps. 



rp_L,3 
rj,_L,2 
rj,_L,l 



^rpx- ~3xy 
^rpz{x^-y' 
^rpx{r^ 
rj^x^ —3xy^ 



rpy 

rpxyzy 

) j^y{r^ 



3yx' 



'5z^)^ 



(BIO) 
(BID 
(B12) 
(B13) 



which take on finite values in the different ordered states. 
These results are summarized in Table |V] As the magnetic 
field breaks time-reversal invariance, the quadrupolar and oc- 
tupolar order parameters may mix linearly in the free energy. 
For example, where is singled out by magnetic field, the 
new terms that enter the free energy are of the form 
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TABLE IV: The symmetry of the different configurations, treating the arrows as polar vectors, or as axial vectors with and without inclusion 
of the time reversal symmetry. In the last column we show the broken symmetry (we assume no magnetic field in the case of the 2:2 collinear 
state and magnetic field along the z direction for 3:1 collinear and for the two canted states). The notation is the same as in the Table lllll with 
the addition of two elements: cr„ is a reflection to a plane perpendicular to the Coo axis (cr^z is also a (t„), while C'2 is a two-fold rotation with 
axis perpendicular to the Coo axis. 

state polar vectors axial vectors spins (axial vectors + time reversal) symmetry broken 

2:2collinear Dooh = Coo®{l,CLa,,J} D^^ Dooh + ODooh 0(3)/(0(2) x 0(1)) = RP^ 

3:1 coUinear Coo.. = Coo x {E, a^} Coo x {£, /} Coo X {E, /} + ecr„Coo X {£, J} 1 

2:2 canted C2« = {£, C2(z), o-a:^, cry-} C^h = {E ,1 ,C'2.{z),a^y) C2h + Qoj:^C2h C00/C2 

3:1 canted Cm = {E, a,,} S2 = {E, 1} S2 + Oa^.Si Coo 
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TABLE V: Classification of rank-three tensor operators according to 
rotational symmetry about a z axis defined by magnetic field. Also 
indicated are the finite values of the order parameters in the 2:2 and 
3:1 canted states. 



order par. 


tensor operators 


2:2 


1:3 









finite 




|yz(3:^-l/^) rpxyz-^ 


finite 


finite 


gi0 


|ya;(r^-5z^) j.y(r-^ -Sz^ ) | 





finite 


1 


jiz(3r^-5z^) 


finite 


finite 



and 

which respect the remaining 0(2) rotational symmetry (more 



precisely, they can mix if the order parameter is finite, ir- 
respectively of the presence of external magnetic field). Mag- 
netic field can therefore strongly modify the symmetry of a 
(primary) multipolar order parameter For a related discus- 
sion, see Ref. ll25n . 



It is not our intention to give a definitive treatment of this 
complex set of coupled order parameters in this paper. How- 
ever we have made a preliminary study of the behavior of 
the rank-three and rank-four tensor order parameters in the 
present model, using the T = theory for a single tetrahe- 
dron embedded in the lattice, and classical Monte Carlo simu- 
lation. We have been unable to identify any higher-order mul- 
tipole which grows faster at a continuous transition than the 
rank-two tensor order parameters given in Section III CI and 
so these retain their tentative assignment as primary order pa- 
rameters. 
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FIG. 23: (Color online) Temperature dependence of (a) nematic or- 
der parameter [Eq. lIlSH. (b) heat capacity [Eq. QUI, (b) and (c) 
the related measure of local correlation Ag"' and its susceptibility 
[Eqs. JlOb and dl IH. and (d) and (e) the global order parameter A|'°''^' 
and its susceptibility [Eqs. ^ and l|9]l] for b = 0.6, h = and a range 
of values of system size. The results are for J3 = —0.02, showing 
both a first order transition into the nematic phase at Tq = 0.518(6) 
(indicated by the vertical dotted line) and a transition into the ordered 
phase at Tjv = 0.42(2) (dashed line). 




FIG. 24: (Color online) Phase diagram for the classical Heisenberg 
antiferromagnet on a pyrochlore lattice, Eq. l[TJ, in applied magnetic 
field h — 0, with additional biquadratic interactions b = 0.6. The 
transition temperature Tn associated with the E-symmetry long- 
range order vanishes as the strength of ferromagnetic third-neighbor 
interactions Jg — s> 0, as determined by Monte Carlo simulation. A 
phase exhibiting nematic order exists above Tn up to a Tq ~ b. 
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FIG. 25: (Color online) (a) Temperature dependence of the primary 
order parameter Q^'^ in the vector-multipole phase. In the limit 
L — > cxD these results extrapolate to a mean field-like behavior 
Q-^'^ ~ VTv - T with Tv ~ 0.091 (solid black points and grey 
line), (b) The secondary order parameter Q^'^ also takes on a finite 
value for T < Ty, but grows more slowly at the transition, (c) Heat 
capacity, showing no measurable singularity at T = Ty- All data are 
for ft. = 5, 6 = 0.1, J3 = 0, with system sizes ranging from L — 4 
to L = 16. 



